Analytic theory of correlation energy and spin polarization in the 2D electron gas 
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Abstract 

We present an analytic theory of the pair distribution function and the ground-state energy in a two-dimensional 
(2D) electron gas with an arbitrary degree of spin polarization. Our approach involves the solution of a zero-energy 
scattering Schrodinger equation with an effective potential which includes a Fermi term from exchange and kinetic 
energy and a Bose-like term from Jastrow-Feenberg correlations. The form of the latter is assessed from an analysis of 
data on a 2D gas of charged bosons. We obtain excellent agreement with data from quantum Monte Carlo studies of 
the 2D electron gas. In particular, our results for the correlation energy show a quantum phase transition occurring 
at coupling strength r s rs 24 from the paramagnetic to the fully spin-polarized fluid. 
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Developments in techniques of molecular beam epitaxy and chemical vapor deposition have allowed since the mid 
1970's the fabrication of semiconductor structures in which the carriers can form a low-density fluid moving in low 
dimensionality [1] . Many of the electron-electron interaction effects become increasingly important as carrier density 
and dimensionality are reduced and the homogeneous electron gas - an assembly of fermions interacting by the e 2 /r 
law and moving in a uniform neutralizing background - provides a primitive model for their study [2]. The model is 
characterized at zero temperature and magnetic field by the carrier density or equivalcntly by the coupling-strength 
parameter r s , which measures the mean intcrparticle distance in units of the Bohr radius as- 

A crucial role in the theory is played by the particle pair distribution function g(r), which determines the ground- 
state energy of the assembly and hence its quantum phase diagram. This function has been the object of several 
numerical studies by quantum Monte Carlo (QMC) methods, starting from the early work of Ceperley and Alder [3] 
on the three-dimensional (3D)clectron gas. In the two-dimensional (2D) case these studies have predicted the critical 
density for Wigner crystallization at r s ss 35 [4,5] and a weakly first-order transition to a fully spin-polarized fluid 
state at r s m 25 [6,7]. 

In a previous work we have presented an analytic theory of g(r) and other ground-state properties of the 3D 
electron gas [8]. The theory was based on a Fermi hypernetted chain approximation (hereafter indicated by the 
FHNC/0 acronym) and yields quantitative agreement with QMC data up to at least r s = 20. The object of the 
present work is to extend the theory of g(r) and the correlation energy to the 2D electron gas. As we shall see, 
much higher sophistication is needed to attain a quantitatively useful theory in lowered dimensionality. We shall 
have to dwell on terms beyond the FHNC/0, which come from low-order elementary diagrams and from three-body 
Jastrow-Feenberg correlations. These effects have been studied theoretically in boson fluids [9,10] and it is natural 
within our theoretical scheme to draw on the 2D gas of charged bosons in dealing with the 2D gas of charged fermions. 
An important role will be played below by the newly available QMC data on the 2D charged-boson gas [11]. 

Our starting point is the energy functional E[n(r)] for an inhomogeneous density profile corresponding to the 
distribution of particles n(r) around an average particle. That is, n(r)= ng(r) , where n = (^rla 2 ^)^ 1 is the average 
particle density, and the corrsponding "external" potential is v(r) = e 2 /r. As usual, the energy functional is the 
sum of an ideal kinetic energy term and of three potential energy terms given by the interaction with the "external" 
potential and by the Hartree and exchange-correlation contributions. As in Ref. [8] we split the ideal kinetic energy 
functional into the sum of the von Weizsacker - Herring term (H 2 /8m) J dr |Vn(r)| 2 /n(r) and a residue Tg[n(r)} , and 
minimize the energy by means of the Hohenberg-Kohn variational principle. The result is a differential equation for 

[-—V 2 + v(r) + W F (r) + W B (r)Wg(r) = . (1) 

This equation is formally exact and describes a zero-energy two-body scattering problem in the many-body fluid, in 
accord with Overhauser's interpretation of g(r) [12]. 

In writing the "scattering potential" in the form used in Eq. (1) we are implicitly thinking of the many-body wave 
function as the product of a Slater determinant times a Jastrow-Feenberg correlation factor. The Fermi potential 
Wf{t) derives from the functional derivative of Ti9[n(r)] and from exchange, and would vanish in a Bose fluid. The 
Bose-like potential Wb(t) contains the effects of correlations (including the long-range Hartree term) and by itself 
would determine g(r) in a Bose fluid. We turn to the latter in order to assess the form of Wb(t) . 

2D fluid of charged bosons. The FHNC/0 and its HNC/0 equivalent for a Bose fluid take 
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where Wsik) is the Fourier transform of Wb(t), defined according to the general expression FT[F(r)] — 
n J drF(r) exp (ik • r). In Eq. (2) S(k) is the structure factor of the fluid, defined so that S(k) = 1 + FT[g(r) - 1]. 
Using Eq. (2) in Eq. (1) with Wf(t) — and self-consistently solving Eq. (1) for g(r), we find quantitative agreement 
with the QMC data of Ref. [11] at r s = 1 and discrepancies already emerging at r s = 5 . 
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Improvements on Eq. (2) can be sought in two directions for a Bose fluid [10,13]. The HNC may be transcended 
by the inclusion of low-order elementary diagrams, and three-body Jastrow-Feenberg correlations may be included. 
The contribution from the first elementary diagrams to the effective Bose potential is 

where 64(h) is given by a fourfold integral (in 2D) over momentum space. The contribution of three-body correlations 
is given by a two-fold integral, 

^W = ^ J rfq5(p)%) U 3(q,P,k){ U (q,p,k) + [ e (p)+ £ (g)] U3 (q,P,k)} (4) 

where p = — (q + k), e(k) = fi 2 k 2 /[2mS(k)] and, with the definition X{k) = 1 — S~ 1 (k), we have w(q, p,k) = 
(H 2 /m) [k • pX(p) + k • qX(q) + p • qX(p)X(q)} and 

(h 2 /2m) [k ■ pX(p)X(fc) + p ■ qX(p)X(q) + k ■ qX(q)X(k)} 
e(fc) + e(p) + £(<?) 

The papers of Smith ei a/. [13] and of Apaja et al. [10] should be consulted for detailed derivations of these equations, 
in which u 3 (q, p, k) is the irreducible three-body vertex. 

The effective scattering potential as self-consistently determined in these approximations for the 2D fluid of charged 
boson at r s = 20 is shown in Fig. 1. In all cases the structure in g(r) which is present in the QMC data of Ref. [11] is 
underestimated at such strong coupling. We have therefore tried to take into account the higher-order terms that are 
missed in these approaches at strong coupling by assuming that they lead to corrections in the scattering potential 
that have roughly the same shape as the low-order terms reported above. In particular, we have found that the choice 

W { B a \r) = WE NC {r) + a(r s )w£\r) (6) 

leads to a satisfactory account of the QMC data on the 2D Bose gas. Here the parameter a is determined by fitting 
the QMC data on the ground-state energy [11] with a relative precision of 10~ 3 . This yields 

a(r s ) = 1 + 5.888 exp (-0.07758 r^- 7923 ). (7) 

The corresponding form of the scattering potential at r s = 20 is also shown in Fig. 1, while Fig. 2 reports a comparison 
of our results for g(r) with the QMC data up to r s = 50. We have clearly achieved fully quantitative agreement over 
this whole range of coupling strength. It may be remarked that nothing is to be gained by further increasing the 
relative precision of the fit of the QMC ground-state energy, possibly because of errors due to the finiteness of the 
QMC sample. 

2D fluid of charged fermions. We carry over Eqs. (4) and (6) in our calculations of the Bose-like term in Eq. (1) 
for the 275 electron gas. With regard to the Fermi term Wf{t) in Eq. (1), we adopt the same criteria in determining 
its form as in Ref. [8]. An important requirement is that Eq. (1) should give the exact fermion fermion distribution 
function when one goes to the weak-coupling limit r s — > 0, when g(r) becomes the Hartree-Fock pair distribution 
function gHF(r)- The Fermi term in the scattering potential is then determined by the Hartree-Fock structure factor 
SHF(k) according to 



h 2 

W F (k) = —FT 
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\f9HF{r) 



[ SH SHF(k) " J [2SHF(k) + ^ W B \k)\s ik)=SHF{k) . (8) 
Here S HF (k) =EA 1 + sgn(a){] Sff F (k)/2 with 

Sff F (k) = (2/tt) [sin- 1 (fc/2fc FCT ) + (k/2k Fa )y/l- {k/2k F M 9(2k F a — k) + 9(k - 2k Fa ) (9) 
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where a = =bl denotes the electron spin, £ = | — nj/n is the spin polarization, fc_p, 



V^l + sgntaX] 1 / 2 /^) 



and is the Heaviside step function. The rationale behind Eq. (8) is as in Ref. [8]: the first term on the RHS ensures 
that the Hartree-Fock limit is correctly embodied into the theory, and the second and third term ensure that the 
Bose-like scattering potential is suppressed for parallel-spin electrons at weak coupling [14]. 

Our numerical results for the pair distribution function of the 2D electron gas in the paramagnetic state (£ = 0) 
and in the fully spin-polarized state (C = 1) are compared with the QMC data of Ref. [15] in Fig. 3. We have clearly 
achieved fully quantitative agreement with the QMC data up to large values of the coupling strength r s . 

We can then confidently calculate the ground-state energy e(r s ,() of the 2D electron gas as a function of r s in 
these two states of magnetization, using an integration over the coupling strength according to the expression 



(in Rydberg units). The results are compared with the QMC data of Attaccalite et al. [7] in Fig. 5. There cleary is 
a discrepancy with the data in the absolute values of the energy, but this affects in essentially the same manner the 
two phases so that we find a transition from the paramagnetic to the fully spin-polarized fluid at r s « 24, in excellent 
agreement with the QMC data. As to the nature of the phase transition, within the accuracy of our calculations it 
could be either a weakly first-order one or a continuous transition occurring in an extremely narrow range of values of 
r s . This is shown in Fig. 4, where we also report our results for the ground-state energy as a function of r s at several 
values of the spin polarization (. 

Conclusions. In summary, we have presented a model which quantitatively predicts the two-body correlations in 
both the 2D charged-boson fluid and the 2D electron gas using as the only input the QMC data for the ground-state 
energy of the boson fluid as a function of the coupling strength r s . The essential physical idea that underlies the model 
is that differences arising from the statistics are disappearing as the fermionic or bosonic fluid is brought into the 
strong coupling regime, where the Coulomb repulsions suppress close encounters of pairs of particles [16]. However, 
exchange between parallel-spin fcrmions must be properly accounted for in the weak-to-intermediate coupling regime. 

The model has allowed us to reproduce the quantum phase transition that has been found to occur in the QMC 
studies of the 2D electron gas, essentially starting from the basic Coulomb Hamiltonian. Within the accuracy of our 
model the transition could be either a weakly discontinuous transition or a continuous one occurring in a very narrow 
range of coupling strength. The latter type of quantum phase transition is found to occur in the 3D electron gas [17]. 
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The effective Bose-like scattering potential (in units of h 2 / \rnr 2 a B )) at 



= 20, as a function of distance r (in 

a B )• The four cases refer to the HNC (Eq. (2)), to the inclusion of a contribution of elementary diagrams 
Wg NC (r) + W B 4 (r) , Eq. (3)) and of three-body correlations (W B (r) = Wg NC (r) + W ( g\r) , Eq. (4)), and 

(6). 
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FIG. 2. The pair distribution function g(r) of the 2D fluid of charged bosons as a function of distance r (in units of r s as) at 
r a = 10, 20, and 50. The theoretical results (full lines) are compared with QMC data by De Palo el al. [11] (dots). The dased 
lines show the results obtained in the HNC/0 scheme. The results have been shifted upwards by 0.4 and 0.8 for clarity. 
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FIG. 3. The pair distribution function g(r) of the 2D electron gas at r s = 1, 5, 10, and 20 in the paramagnetic state (top 
panel) and in the fully spin-polarized state (bottom panel), as a function of distance r (in units of r s a B )• The theoretical 
results (full lines) are compared with QMC data by Gori-Gorgi et al. [15] (dots). The dashed lines show the results obtained in 
the FHNC/0 scheme. The results have been shifted upwards by 0.4, 0.8 and 1.2 for clarity. 
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FIG. 4. The ground-state energy of the 2D electron gas (in Rydberg units, referred to the Madelung energy em = — 2.2122/r s 
and multiplied by r z J 2 ) as a function of the coupling strength r s . The full lines show the theoretical results for the paramagnetic 
state and the fully spin-polarized state, while the dots report QMC data from Attaccalite el al. [7]. The other three curves are 
theoretical results for states of partial spin polarization, as indicated in the figure. 
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